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The massive end of the luminosity and stellar mass 
functions: Dependence on the fit to the light profile 
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ABSTRACT 

In addition to the large systematic differences arising from assmnptions about the stel- 
lar mass-to-light ratio, the massive end of the stellar mass function is rather sensitive 
to how one fits the light profiles of the most luminous galaxies. We quantify this by 
comparing the luminosity and stellar mass functions based on SDSS cmodel magni- 
tudes, and PyMorph single-Sersic and Sersic-Exponential fits to the surface brightness 
profiles of galaxies in the SDSS. The PyMorph fits return more light, so that the pre- 
dicted masses are larger than when cmodel magnitudes are used. As a result, the total 
stellar mass density at z ~ 0.1 is about 1.2 x larger than in our previous analysis of the 
SDSS. The differences are most pronounced at the massive end, where the measured 
number density of objects having Af* > 6 x lO^'^M© is ■^ 5x larger. Alternatively, at 
number densities of lO^^Mpc^'^, the limiting stellar mass is 2x larger. The differences 
with respect to fits by other authors, typically based on Petrosian-like magnitudes, 
are even more dramatic, although these differences are sometimes masked by large 
differences in the assumed M^jL (even after scaling to the same IMF). Our results 
impact studies of the growth and assembly of stellar mass in galaxies, and of the re- 
lation between stellar and halo mass, so we provide simple analytic fits to these new 
luminosity and stellar mass functions and quantify how they depend on morphology, 
as well as the binned counts in electronic format. 
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1 INTRODUCTION 

The brightest, most massive galaxies have been the object of 
much study. Recent work has emphasized the importance of 
using a good parametrization of the abundance at the bright, 
massive end if one is interested in using Halo Model based 
abundance matching techniques, or extreme value statistics, 
to understand their origin (e.g. Paranjape & Sheth 2012). 
A few years ago Bernardi et al. (2010) noted that the most 
luminous galaxies were more abundant than expected from 
the most commonly used parametrizations of the luminos- 
ity function. They also pointed out that, when converted to 
a stellar mass function, this mis-match was important for 
models which use the observed abundance and its evolution 
to constrain the issue of whether these objects were assem- 
bled via major or minor mergers. However, they also showed 
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that the conversion from <i>(L') to (^(M-t) is rather sensitive 
to the assumed stellar mass-to-light ratio, for which, as we 
show below, there is still no consensus. 

Bernardi et al. (2010) used luminosities estimated from 
the cmodel magnitudes of the Sloan Digital Sky Survey 
(hereafter SDSS). Although for the brightest galaxies they 
tended to return more light than the more commonly used 
estimates based on the Petrosian radius, they are still a 
poor-man's best guesstimate for the total light if the sur- 
face brightness distribution of the objects follows neither a 
pure exponential disk nor a deVaucouleur's profile (Bernardi 
et al. 2007). Recently, Meert et al. (2013a, b) have performed 
more careful Sersic-bulge -I- exponential disk decompositions 
of these objects. These typically return even more light than 
the cmodel magnitudes (e.g. Bernardi et al. 2013), so the 
main purpose of the present note is to show how this im- 
pacts estimates of the luminosity and stellar mass functions 
at the bright end. As one might expect, the effect is dramatic 
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- at least as dramatic as the choice of Mt/L. Therefore, a 
related goal of the present work is to separate out the effect 
on 4>{M*) of how the luminosity was estimated from that of 
M,/L. 

The next Section describes our sample, shows the lu- 
minosity and stellar mass functions, quantifies how they 
depend on the fit to the light profile and provides simple 
fitting formulae which quantify our results as well as the 
binned counts in electronic format. For reasons described 
in Bernardi et al (2010), all our M,/L estimates assume 
a Chabrier (2003) IMF. In Section 3 we show that, even 
though a number of recent works have made this same choice 
for the IMF (Baldry et al. 2012; Moustakas et al. 2013), 
they still have Af,/L values which are very difi'erent from 
ours (i.e., Bernardi et al. 2010), from one another, and from 
earlier work (Bell et al. 2003). While our results allow one 
to easily account for the dependence on the light profile 
(e.g. use Sersic instead of Petrosian magnitudes), the ques- 
tion of which Mi,/L estimate is most appropriate is beyond 
the scope of this work, and deserves further study. In Sec- 
tion 4 we show how the luminosity and stellar mass func- 
tions depend on morphological type, where the type is deter- 
mined by the Bayesian Automated Classification scheme of 
Huertas-Company et al. (2011). A final section summarizes. 

When converting from apparent brightnesses to lumi- 
nosities, we assume a spatially fiat background cosmology 
dominated by a cosmological constant, with parameters 
{Q.m, ^a) = (0.3, 0.7), and a Hubble constant at the present 
time of Ho = 70 km s~^Mpc"\ 



2 LUMINOSITY AND STELLAR MASS 
FUNCTIONS 

2.1 The sample 

To provide a direct comparison with previous work, we have 
selected the same sample as Bernardi et al. (2010). I.e., 
about 260,000 SDSS galaxies having 14.5 < m,.pot < 17.7. 
We ran PyMorph (Vikram et al. 2010, Meert et al. 2013a) on 
these objects. This is an algorithm which uses GALFIT (Peng 
et al. 2002) to fit seeing convolved 2- dimensional Sersic -I- 
exponential models to the observed surface brightness pro- 
files of galaxy images. Results from extensive tests indicate 
that the algorithm is accurate (Meert et al. 2013a). PyMorph 
sometimes fails to converge to an answer; this happens about 
2% of the time, but because this fraction is independent of 
magnitude, it does not affect our completeness, other than 
by a small overall scaling. We then computed k- and evo- 
lution corrections for each object following Bernardi et al. 
(2010), and hence, luminosities. 



2.2 Dependence on assumed surface brightness 
profile 

The magnitudes and half-light radii output by PyMorph de- 
pend on the model which is fit. E.g., fitting what is really 
a two-component image with a single deVaucouleurs profile 
will generally underestimate the total light. On the other 
hand, the total light associated with the best-fit single Sersic 
or a two-component Sersic-bulge -I- exponential-disk model, 
is less-biased from its true value (e.g. Bernardi et al. 2007; 
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Figure 1. Difference between PyMorph Sersic fits and SDSS Pet- 
rosian, SDSS cmodel, PyMorph SerExp, and Sersic fits from Simard 
et al. (2011) (bottom to top), for galaxies in the sample selected 
by Bernardi et al. (2010). Petrosian magnitudes are always the 
faintest, whereas single Sersic-based magnitudes tend to be the 
brightest. Dotted lines around PyMorph (Ser)-Simard (Ser) show 
the 16th and 84th percentiles of the distribution; these are similar 
to the scatter around the median for the other curves. 



Bernardi et al. 2013). Meert et al. (2013a) and Bernardi et 
al. (2013) have also shown that in objects brighter than L, , 
fitting a two-component Sersic + exponential model to what 
is really just a single Sersic results in a noisier recovery of 
the input parameters, but these are not biased. On the other 
hand, fitting a single Sersic to what is truly a two-component 
system results in significant biases. 

Although the Sersic -I- exponential model is more accu- 
rate, the Sersic fit is often performed on real data when it 
is believed that the resolution and S/N are such that it is 
unlikely to recover a robust two-component fit. Therefore, 
since either of these models are expected to be more realistic 
than a single deVaucouleurs model, we will use both in what 
follows. 

The estimates of the total light from a Sersic or Sersic -f- 
exponential model are generally larger than those based on 
the cmodel magnitudes output by the SDSS pipelines (e.g. 
Bernardi et al. 2007; Hill et al. 2011; Bernardi et al. 2013; 
also see Mosleh, Williams & Franx 2013), and both are larger 
than the SDSS Petrosian magnitudes. Figure 1 illustrates 
that this difference can be large. (The cmodel magnitudes 
were corrected for the SDSS sky subtraction problems as 
described in Bernardi et al. 2010, but this is a small correc- 
tion.) Although the strong dependence on the assumed light 
profile is what has motivated our study, it is reasonable to 
ask if these differences are indeed larger than those associ- 
ated with different pipelines which fit the same model. We 
address this in the next subsection. 



2.3 Dependence on pipeline 

As a check, we have also used luminosities from the single- 
Sersic-based photometric reductions of Simard et al. (2011). 
Figure 1 shows that these are in good agreement with 
PyMorph except for a small offset (~ 0.05 mags), although 
the differences become large at the bright end. See Fig- 
ures Al and A2 in Bernardi et al. (2013) and discussion on 
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Figure 2. Comparison of single Sersic reductions for tlie SDSS galaxies in common to Simard et al. (2011), Kelvin et al. (2012) and 
PyMorph (our notation Mg^^-iORe reflects tlie fact that Kelvin et al. truncate the profile at lOiJe)- Red triangles show a similar analysis 
if one restricts to objects with z < 0.06 as done in Baldry et al. (2012); this shallower volume does not probe the highest luminosities 
that are of most interest here. Dotted lines show the 16th and 84th percentiles of the distribution. 



sky estimates in Meert et al. (2013a,b) for why we believe our 
estimates are less biased. In any case, these differences are 
small compared to PyHorph-cmodel or PyMorph-Petrosian. 

The PyMorph and Simard et al. luminosities come from 
integrating the fitted profile to infinity. Other authors trun- 
cate, typically at some multiple of the half-light radius. 
For example, the analysis of galaxies in the GAMA survey 
Kelvin et al. (2012) truncates the fits at lORe- Using the 
GAMA DRl data release (Driver et al. 2011), we compare 
PyMorph, Simard, and GAMA values for the 7335 galaxies 
for which all three reductions are available. (This sample 
is set by the fact that the GAMA DRl covers 100 sq.deg. 
of the SDSS. GAMA has 10750 matches with the DR7 
SDSS spectroscopic galaxy sample, of which 7335 galax- 
ies are in the Bernard! et al. 2010 sample we study here.) 
Figure 2 compares PyMorph, Simard, and GAMA values for 
the single-Sersic magnitude. The bottom right panel shows 
that the truncation matters at the level of 0.05 mags only 
at Mr < —22. But otherwise, if truncated similarly, then 
GAMA and PyMorph are in good agreement at Mr > —22, 
whereas PyMorph returns significantly more light than the 
other two at the bright end. We believe the differences at 
the bright end are similar in origin (i.e. sky subtraction is- 
sues) to those with respect to Simard et al. (see Meert et al. 
2013a,b; Bernardi et al. 2013). 



2.4 The luminosity function 

For each of the estimates of the total light shown in Figure 1, 
we estimated the luminosity function as in Bernardi et al. 
(2010) using the Knax method of Schmidt (1968). (I.e., we 
weighted each galaxy using l/ymax(ipct), where Knax is the 



maximum comoving volume within which the object could 
have been included in the sample, accounting for both the 
bright and faint magnitude limits.) 

Figure 3 shows the luminosity functions for the SDSS 
Petrosian and cmodel magnitudes (corrected for the SDSS 
sky subtraction problems as described in Bernardi et al. 
2010), and SerExp and Sersic magnitudes (from PyMorph). 
Although the difference between the Petrosian and cmodel 
magnitudes has been known for some time, the fact that 
single-Sersic based counts lie substantially above those based 
on the SDSS outputs has only recently begun to attract at- 
tention. For example, the GAMA based results of Hill et 
al. (2011) point to this difference, but because GAMA cov- 
ers a substantially smaller volume than the SDSS, it does 
not probe the high luminosity end which is of most interest 
here. Our PyMorph reductions, which are in good agreement 
with Hill et al. at Mr > —23, show that at Mr < —23 the 
difference with respect to cmodel counts is dramatic indeed. 

The PyMorph-based counts are in good agreement with 
those which use the Simard et al. (2011) single-Sersic re- 
ductions, except at luminosities brighter than Mr ~ —24, 
where PyMorph tends to be brighter (c.f. Figure 1), so the 
PyMorph luminosity function shows more high luminosity ob- 
jects. This agreement illustrates that our finding that single- 
Sersic fits return substantially more objects in the high lumi- 
nosity tail than do cmodel magnitudes is robust to changes 
in the reduction pipeline. 

The solid curves show the result of fitting 
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with X = Lto the counts. The associated luminosity density 
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Figure 3. SDSS Main galaxy luminosity function based on 
Petrosian, cmodel, single Sersic from Simard et al. (2011) and 
PyMorph SerExp and Sersic magnitudes (bottom to top at Mr = 
—24). Smooth curves show the result of fitting equation (1) to 
the counts; associated best-fit parameter values are given in Ta- 
ble 1. For the Petrosian and cmodel magnitudes, the curve shown 
is that reported by Bernardi et al. (2010) on the basis of fitting 
to Mr < —20. The Petrosian and Sersic based fits of Bell et al. 
(2003) and Hill et al. (2011), respectively, are also shown for com- 
parison. 
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Figure 4. Same as previous figure, but now for the associated 
stellar mass functions. Recent stellar mass functions from Baldry 
et al. (2011; based on Sersic magnitudes) and Moustakas et al. 
(2013; based on cmodel magnitudes) are also shown. All stellar 
masses assume a Chabrier IMF. 



ispx = (t>cX,r[{l + a)/l3]/r[a/l3] + (l,-,X^r[l + 'y]. The first 
term in equation (1) is the same functional form as that used 
by Bernardi et al. (2010); the second is required to fit the 
slight bump at the faint-end. The parameters which yield the 
best-fit are given in Table 1. Note that the value of X, is 
not as intuitive as is its mean value Xt T[{l + a)/l3]/T[a/l3]. 
The observed distributions shown here have been broad- 
ened slightly by measurement errors. Bernardi et al. (2010) 
showed how to modify the analog of equation (1) so as to es- 
timate the parameters of the intrinsic distribution, but that, 
in practice, the difference between the intrinsic and observed 
broadened distributions is small - much smaller than the dif- 
ference between the PyMorph and cmodel counts, so here we 
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Figure 5. Similar to previous figure, but now showing cumula- 
tive rather than differential counts. Top and bottom panels show 
number and stellar mass density respectively. To facilitate com- 
parison with previous work we show the fit of Bell et al. (2003). 



show the results of the observed distribution not the intrinsic 
one. 



2.5 The stellar mass function 

Figure 4 shows the associated stellar mass functions. In all 
cases, A'/« was estimated from the luminosity and the cmodel 
g — r color assuming the Chabrier (2003) IMF as described 
in Bernardi et al. (2010). We use the cmodel color because 
the main goal of this paper is to study the effect on 0(M») 
from changes in L. By using cmodel colors, we are ensur- 
ing that our M^/L estimates for each object are the same 
as in Bernardi et al. (2010); however, the L estimate for 
each object differs (Petrosian 7^ cmodel 7^ PyMorph). Notice 
again that the PyMorph-based estimates (as well as those 
from Simard et al.) lie well-above the Petrosian and cmodel 
ones. (Of course, if the stellar population models used to es- 
timate MtjL are incorrect, or if the IMF is mass-dependent, 
then this will modify the results. See Section 3 for compari- 
son with other work.) 

The estimate from Baldry et al. (2011) lies below all 
the others. This is remarkable because it is based on the 
GAMA-Sersic reductions, and we have already seen that the 
associated <^(I/) (from Hill et al. 2011) is in good agreement 
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Table 1. Parameters of (j>{Lr) (top rows) and (f>(M-^) (bottom rows) derived from fitting equations (1) to the observed counts based on 
different magnitudes. 
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with that based on PyMorph. Therefore, the difference in 
<j){Mt) must be entirely due to Mt/L, even though Baldry 
et al. also assume a Chabrier IMF. We discuss this more in 
Section 3. 

We think it is interesting to present our results in a 
format which highhghts just how much the PyMorph-bascd 
values differ from other work (we use Bell et al. 2003 for com- 
parison). Figure 5 shows cumulative (rather than differen- 
tial) counts, both for number and stellar mass- weighted den- 
sity. The number counts at the mass scale above which the 
number density of objects is lO^^Mpc^"^ is larger by a factor 
of ~ 2 compared to the cmodel-based counts (a factor of ~ 3 
compared to Bell et al.). Alternatively, at M, = 6 x 10^^ M©, 
the PyMorph counts lie a factor of ~ 8 above those based on 
cmodel magnitudes (a much larger factor above Bell et al.). 
For the mass-weighted counts the corresponding discrepan- 
cies at lO^MoMpc"^ or 6 X lO^M© are similar or slightly 
larger. 

To make our results simple to use, in addition to Ta- 
ble 1, which reports the parameter values associated with 
the best-fits, we have made the binned counts available in 
the electronic versions of Tables 2 and 3. 



3 COMPARISON WITH PREVIOUS WORK 

Recently, it has become fashionable to concentrate more on 
(j>{Mt) than 4>{L). Unfortunately, this combines two very 
different types of uncertainty: that associated with the to- 
tal light, and the other associated with the stellar mass-to- 
light ratio Mt/L (e.g. mismatched stellar templates, mass- 
dependence of the IMF, etc.). To illustrate this, we compare 
our results with two of the most recent determinations of 
(j){Mt): those of Baldry et al. (2011) and Moustakas et al. 
(2013). Although we all assume a Chabrier IMF, their esti- 
mates of (jf>(M,) are more similar to one another than they 
are to ours. However, as we argue below, this implies large 
differences in their M*/L values since they used different 
luminosities to get M, . 

The Baldry et al. analysis is based on single Sersic fits to 
the light profiles from galaxies with z < 0.06 in the GAMA 



survey reported by Kelvin et al. (2012). Of the 7335 GAMA 
matches in the Bernardi et al. (2010) sample which we are 
studying here, only 1612 have z < 0.06. The red triangles in 
Figure 2 show that, for these 1612 galaxies, the single Sersic 
PyMorph estimates of the total light are in good agreement 
with those derived by Kelvin et al., and used by Baldry et 
al. However, the redshift cut eliminates most of the high 
luminosity objects which are of most interest to our paper. 

We have checked explicitly and found that the GAMA- 
based luminosity function for these z < 0.06 objects is ac- 
tually in good agreement with that based on PyMorph single 
Sersic reductions when restricted to z < 0.06, and this is 
in good agreement with that for the full GAMA sample of 
Hill et al. (2011) shown in Figure 3. (Of course, the smaller 
volume means the error bars are large, and the compari- 
son is effectively limited to abundances greater than about 
10~* Mpc~^ dex~^.) However, our PyMorph-based 0(M,) es- 
timate lies well above that of Baldry ct al., so we conclude 
that our M,/L values must be larger than theirs. The Baldry 
et al. estimate is much closer to, though slightly below, the 
estimate of Bell et al. (2003, scaled to a Chabrier IMF). 
However, since the Baldry et al. Sersic-based <j}{L) estimate 
agrees with ours, and this last lies well above the 0(1/) as- 
sociated with the Petrosian magnitudes which were used by 
Bell et al. (see our Figure 3) we conclude that the Baldry et 
al. Mt/L values nmst also be smaller than those of Bell et 
al. 

Moustakas et al. use cmodel magnitudes. Since their 
sample is essentially the same SDSS sample as ours, we ex- 
pect their luminosity function to agree with that of Bernardi 
et al. (2010). However, Figure 4 shows that although their 
stellar mass function is similar to our cmodel-derived 0(Mt) 
at low masses, it is different at the high end, indicating that 
the Moustakas et al. Mt/L ratios are smaller than ours. (We 
suspect that this must be related to the choice of template 
used to estimate Mt/L at higher masses, for reasons given 
in Fig. 22 and associated discussion of Bernardi et al. 2010.) 
On the other hand, despite differences at lower masses, the 
Moustakas et al. (j>{Mt) is reasonably well approximated by 
the Bell et al. estimate at higher masses. Since the cmodel 
luminosity function is different from the Petrosian one at 
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these high masses, we conclude that the Moustakas et ah 
Mt/L values must be smaller than those of Bell et al., and 
different again from those of Baldry et al. (who used Sersic 
rather than cmodel magnitudes). 
Therefore, 

• by comparing the PyMorph Sersic-based (/)(M«) with 
that based on Simard et al. our Figure 4 quantifies how 
differences due to the uncertainties in a given light profile 
fit (mainly due to sky subtraction issues, see Meert et al. 
2013b) affect </>(M.); 

• by comparing our Petrosian-based (j){A4f) to that 
based on cmodel, SerExp and Sersic fits, our Figure 4 quan- 
tifies the effect of fitting different light profile models; 

• by comparing our Petrosian based <j>{Mt) to the Bell 
et al. fit, our cmodel-based estimate to that from Moustakas 
et al., or our Sersic (or Simard et al.) based 0(M*) to Baldry 
et al, our Figure 4 quantifies how systematic differences in 
M*/L affect (/)(A'/.). 

This shows that the effects on (f>{M,) of using the total lu- 
minosity computed from different fits to the light profile 
are dramatic; it is important to specify how the light profile 
was fit when reporting a luminosity or stellar mass function. 
That said, a number of authors now agree that Petrosian and 
cmodel magnitudes underestimate the total light (Simard et 
al. 2011; HiU et al. 2011; Kelvin et al. 2012; Bernardi et al. 
2013; Meert et al. 2013b; Mosley et al. 2013). The differences 
between these and our PyMorph Sersic-based magnitudes are 
larger than those between the PyMorph and Simard et al. 
Sersic magnitudes (see Figure 1). So we believe that future 
estimates of (j)(M*) should be based on Sersic or SerExp 
rather than cmodel magnitudes, and, especially at the high 
mass end, Petrosian magnitudes should not be used. 



4 DEPENDENCE ON MORPHOLOGY 

We have combined our PyMorph SerExp reductions with the 
Bayesian Automated morphological Classifier of Huertas- 
Company et al. (2011) to determine how the luminosity 
and stellar mass functions depend on morphology. This al- 
gorithm returns the probability that an object is one of four 
types (for each object, the sum of the four probabilities is 
unity). Therefore, to estimate the luminosity function, we 
simply weighted object j by pj (type) /Vma.x{L^^^) . 

In practice, there are a number of faint objects for which 
Pi(E) < 0.15. These can dominate over the counts from simi- 
larly faint objects for which Pj(E) > 0.85. If these low values 
of p axe simply the result of errors in the BAG algorithm, 
then these will wrongly boost the luminosity function at the 
faint end. To check the magnitude of the effect, we have set 
to zero all values of p < 0.15, and reassigned the weight to 
the types which remain with contribution proportional to 
the nonzero remaining values of p such that the sum over 
the four p values (some of which are now zero) is still unity. 
This reduces the counts of faint ellipticals and luminous spi- 
rals by a factor of about two: it is these modified counts 
which we show in Figure 6. 

The binned counts are given in Tables 2 and 3, which 
are provided in convenient electronic format in the online 
version of the journal. (Note that summing up the luminos- 
ity functions of the different morphological types does not 
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Figure 6. Morphological dependence of the luminosity (top) and 
stellar mass (bottom) functions. In all cases, the luminosities are 
based on Sersic magnitudes, stellar masses assume a Chabrier 
IMF, and morphologies are based on the BAG method of Huertas- 
Company et al. (2011). 



quite give the luminosity function of the full sample. The 
small discrepancy arises because we were unable to find BAG 
classifications for a small fraction of the objects in our sam- 
ple.) Figure 6 confirms the well-known trend for early-types 
(E and SO) to dominate the high-mass end, and later-types 
(Sa to Sd) to dominate at lower masses. 

We have compared these estimates with those of Nair & 
Abraham (2010), who used the T-Type classification using 
the modified RG3 classifiers. We set E (T= -5 and T= -4), 
SO (T= -3, T= -2 and T= -1), Sa (T= 0, T= 1 and 
T= 2), Sb (T= 3 and T= 4), and Scd (T= 5, T= 6 and 
T= 7). Our E and SO counts are in quite good agreement, 
though, at the bright end, they tend to classify objects as 
E rather than SO when BAG divides its weights between E 
and SO. Differences are slightly more pronounced for later 
types: our Scd counts are in good agreement at the faint end, 
but lie substantially above theirs at high luminosities (where 
the counts are falling exponentially). On the other hand, our 
Sab counts lie about a factor of two below theirs at high L, 
suggesting that BAG assigns some weight to the extreme 
Scd classification when Nair et al. choose the extremes. 
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Table 2. The binned (f>(Mr) counts for the full sample and when weighted by the probability of a given morphological type. X = 
Mr [mag] and Y = Logio</'(Mr) [Mpc~^ dex~^]. Four electronic tables are provided in this format based on the type of magnitudes; 
PyMorph Sersic (LF-Ser.dat), PyMorph ScrExp (LF-SerExp.dat), Slmard Sersic (LF-Ser-Simard.dat) and cmodcl from Bernardi et al. 
(2010) (LF-cmodel.dat). 



X 


Y (All) 


Y wP(Ell) 


Y wP(SO) 


Y wP(Sab) 


Y wP(Scd) 


-17.700 


-2.350±0.065 


-4.030±0.065 


-3.209±0.065 


-2.708±0.065 


-2.706±0.065 



Table 3. The stellar mass function 0(M*) of the full sample, and when weighted by the probability of being a given morphological 
type. X = LogioM* [Mq] and Y=Logio[(ln 10) M, </>(M*)] [Mpc~^ dex"-']. Four electronic tables arc provided in this format based 
on the type of magnitudes: PyMorph Sersic (MsF-Ser.dat), PyMorph SerExp (MsF-SerExp.dat), Simard Sersic (MsF-Ser-Simard.dat) and 
cmodel from Bernardi ct al. (2010) (MsF-cmodcl.dat). 



X 



Y (All) 



Y wP(Ell) 



Y wP(SO) Y wP(Sab) Y wP(Scd) 



9.050 



-2.012±0.053 -3.884±0.053 -2.933±0.053 -2.376±0.053 -2.339±0.053 



5 DISCUSSION 

PyMorph Sersic or Sersic+exponential based estimates of the 
total light of a galaxy are larger than those based on SDSS 
pipeline outputs (Petrosian or cmodel; Figure 1). As a result 
our PyMorph-based luminosity and stellar mass functions are 
rather different from previous work: they have more light 
and mass at the bright, massive end (Figures 3 and 4). Our 
Sersic-based luminosity functions, from both PyMorph and 
Simard et al. (2011) reductions, are very different from the 
Sersic-based analysis in Blanton et al. (2003) who found no 
difference between the Petrosian and Sersic-based luminos- 
ity functions. However, the Sersic parameters used by Blan- 
ton et al. were estimated from a 1-dimensional radial surface 
brightness profile, measured in ~ 5 — 10 azimuthally aver- 
aged annuli. This procedure is expected to be significantly 
less accurate than the 2-dimensional fits to the whole galaxy 
image performed by PyMorph and Simard et al. Therefore, 
we believe that these more recent Sersic reductions (as well 
as those of Kelvin et al. 2012 for the subset of objects in the 
GAM A survey) supercede those of Blanton et al., making 
obsolete the Blanton et al. conclusions about the Sersic lu- 
minosity function being similar to the Petrosian-based one. 

Our fits indicate that the luminosity density at z ~ 0.1 
is about 10% larger than our previous work with cmodel 
magnitudes (Table 1), and about a factor of two larger than 
when based on Petrosian magnitudes (Figure 5). This differ- 
ence is driven by the most luminous objects which are pre- 
dominantly quiescent early- type galaxies (Figure 6). Since 
a number of authors now agree that SDSS pipeline mag- 
nitudes underestimate the true luminosity, and these more 
recent algorithms are in reasonably good agreement with 
the PyMorph reductions used by us (Figure 2), we conclude 
that there is now good agreement that the bright end of 
the luminosity function may be substantially brighter than 
Petrosian magnitudes suggest. 

As one might expect given our analysis of the lumi- 
nosity function, our (Chabrier IMF-based) stellar mass den- 
sities at z ~ 0.1 are about 25% larger than our previous 
work (Bernardi et al. 2010) with cmodel magnitudes (Ta- 



ble 1), which was itself considerably larger than when based 
on the Petrosian magnitudes which are often used for this 
purpose (Figure 5). As a result, our estimates have impli- 
cations for studies of the evolution of the star formation 
rate, the growth of stellar mass in galaxies, the processes by 
which this mass was assembled, and Halo Model analyses of 
the Mstar — Mhaio relation. For example: 

• Our higher stellar mass density at z ~ resolves 
the tension with respect to the total mass density inferred 
from the integrated star formation rate (SFR), as noted in 
Bernardi et al. (2010). 

• A higher number density of massive galaxies in the local 
Universe allows for a higher incidence of major (in addition 
to minor) mergers in driving the stellar mass growth of the 
most massive central galaxies at late times (e.g., Bernardi 
et al. 2011a,b; Shankar et al. 2013). This conclusion rests, of 
course, on the quality of the determination of the high red- 
shift stellar mass function, a task which we expect to be even 
more challenging than in the local Universe. In this respect, 
it is interesting that the z ^ 1 counts at M* > lO^^'^M© 
of CaroUo et al. (2013) are in rather good agreement with 
our 2: ~ (Sersic-based) estimate, strongly limiting merger 
rates at the high mass end. 

• A higher number density at high masses would better 
match the stellar and dynamical mass functions, possibly re- 
ducing the need for a strong mass-dependent variation of the 
IMF (see, e.g.. Fig. 23 in Bernardi et al. 2010), although the 
true extent of the latter statement relies on accurate dynam- 
ical mass measurements with appropriate effective radiuses 
and structure constants. We plan to address this separately. 

• The stellar mass function in the local Universe is one of 
the fundamental ingredients in popular semi-empirical mod- 
els for populating haloes with galaxies, such as the halo oc- 
cupation and abundance matching techniques (e.g., Cooray 
& Sheth 2002; Berlind & Weinberg 2002; Vale & Ostriker 
2004; Shankar et al. 2006; Zehavi et al. 2011; Leauthaud et 
al. 2012; Moster et al. 2012). If the more massive galaxies 
are more abundant, then Halo Model analyses will assign 
them to lower mass halos. Since lower mass halos are less 
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strongly clustered, we expect the most massive galaxies to 
be less strongly clustered than current models assume. 

• This would also imply that the median baryon fraction 
at the high mass end may be significantly higher than pre- 
viously thought. This may pose serious questions about the 
impact of feedback from active galactic nuclei, in the quasar 
and/or radio modes (e.g., Granato et al. 2004; Croton et 
al. 2006; Silk & Mamon 2012); even more so, if one consid- 
ers that the baryon fraction of ellipticals at their formation 
epoch [z > 1) must have been even higher. 

• More stellar mass at the high mass end directly impacts 
studies of how the stellar fraction compares with that in the 
gas detected by X-ray and Sunyaev-Zeldovich experiments. 
E.g., by decreasing the halo mass one should associate with 
a given stellar mass, it potentially reduces the discrepancy 
shown in Figure 9 of Ade et al. (2013). 

• These effects on the stellar to halo mass mapping go in 
the same direction as those from plausible changes to the 
IMF at the massive end (Bernardi et al. 2010 and references 
therein), so the data may be indicating that a major revision 
of the results on the galaxy-halo mapping at the high mass 
end, is called- for. 

Since a number of groups are currently engaged in such stud- 
ies, we have provided the luminosity and stellar mass func- 
tions shown in Figures 3, 4 and 6 in a convenient electronic 
tabular format in the online version of the journal (see Ta- 
bles 2 and 3). 

We caution that our estimates of the amount of stel- 
lar mass in the most massive objects are at least 2x larger 
than other recent determinations (e.g. Baldry et al. 2012 
and Moustakas et al. 2013), even though we argued that 
our estimates of the luminosity function are in much bet- 
ter agreement than this (;/>(M»)-based number suggests. Sec- 
tion 3 shows that their M* values are similar only because 
they assume very different L and M^/L values from one 
another. Therefore, until the field converges on what it be- 
lieves to be reliable M-^/L estimates for the highest mass 
objects (see Mitchell et al. 2013 for ongoing discussion of 
this point), we believe our results argue against calibrating 
models to published stellar mass functions: calibrating to the 
luminosity function is more robust, as it does not combine 
differences in L and M,/L into a single number M*. 

When using our results, it is important to bear in mind 
that the PyMorph estimates assume that the galaxy is either 
a single Sersic, or a two-component Sersic+Exponential sys- 
tem. These both allow for substantial light beyond the core 
of the image, so there is some question as to just what it is 
that the profiles are fitting. The most massive objects tend 
to be BCGs, so some of the excess light returned by Sersic 
and/or SerExp fits is likely to contain intracluster light. If 
so, then our estimates of the total stellar mass are appro- 
priate for models which associate the ICL with the central 
galaxy. This may even be physically reasonable, since most 
of the accretion and stripping which occured as the clus- 
ter assembled likely happened during accretion onto what is 
now the central object. 
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